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1. INTRODUCTION 



Though commonly umecognized, a superconducting BCS condensate consists of 
equal numbers of two-electron (2e) and two-hole (2h) Cooper pairs (CPs). A com- 
plete boson-fermion (statistical) model (CBFM), however, is able to depart from this 
perfect 2e-/2h-CP symmetry and yields [1] robustly higher T^s without abandon- 
ing electron-phonon dynamics mimicked by the BCS/Cooper model interaction Vk,k' 
which is a nonzero negative constant —V, if and only if single-particle energies e^, ek' 
are within an interval [max{0, ^ — huj d} , IJ' + d] where fi is the electron chemi- 
cal potential and is the Debye frequency. The CBFM is "complete" only in the 
sense that 2h-CPs are not ignored, and reduces to all the known statistical theo- 
ries of superconductors (SCs), including the BCS-Bose "crossover" picture but goes 
considerably beyond it. 

Boson-fermion (BF) models of SCs as a Bose-Einstein condensation (BEC) go 
back to the mid-1950's [2-5], pre-dating even the BCS-Bogoliubov theory [6-8]. Al- 
though BCS theory only contemplates the presence of "Cooper correlations" of single- 
particle states, BF models [2-5, 9-17] posit the existence of actual bosonic CPs. In- 
deed, CPs appear to be universally accepted as the single most important ingredient 
of SCs, whether conventional or "exotic" and whether of low- or high-transition- 
temperatures Tc. In spite of their centrality, however, they are poorly understood. 
The fundamental drawback of early [2-5] BF models, which took 2e bosons as analo- 
gous to diatomic molecules in a classical atom-molecule gas mixture, is the notorious 
absence of an electron energy gap A(T). "Gapless" models cannot describe the su- 
perconducting state at all, although they are useful [16,17] in locating transition 



temperatures if approached from above, i.e., T > Tc- Even so, we are not aware of 
any calculations with the early BF models attempting to reproduce any empirical 
values. The gap first began to appear in later BF models [9-14]. With two [12, 13] 
exceptions, however, all BF models neglect the effect of hole CPs accounted for on an 
equal footing with electron CPs, except the CBFM which consists of both bosonic CP 
species coexisting with unpaired electrons, in a ternary gas mixture. Unfortunately, 
no experiment has yet been performed, to our knowledge, that distinguishes between 
electron and hole CPs. 

The "ordinary" CP problem [18] for two distinct interfermion interactions (the 
5- well [19, 20] or the Cooper/ BCS model [6, 18] interactions) neglects the effect of 2h 
CPs treated on an equal footing with 2e [or, in general, two-particle (2p)] CPs. On 
the other hand. Green's functions [21] can naturally deal with hole propagation and 
thus treat both 2e- and 2h-CPs [22, 23]. In addition to the generalized CP problem, a 
crucial result [12, 13] is that the BCS condensate consists of equal numbers of 2p and 
2h CPs. This was already evident, though widely ignored, from the perfect symmetry 
about e = of the well-known Bogoliubov [24] v^{e) and u^{e) coefficients, where e 
is the electron energy. 

Here we show: a) how the crossover picture TcS, defined self-consistently by 
both the gap and fermion-number equations, requires unphysically large couplings 
(at least for the Cooper/BCS model interaction in SCs) to differ significantly from 
the Tc from ordinary BCS theory defined without the number equation since here 
the chemical potential is assumed equal to the Fermi energy; how although ignoring 
either 2h- or 2e-CPs in the CBFM b) one obtains the precise BCS gap equation for 
all temperatures T, but c) only half the T = BCS condensation energy emerges. 
The gap equation gives A(T) as a function of coupling, from which Tc is found as the 
solution of A(Tc) = 0. The condensation energy is simply related to the ground-state 
energy of the many-fermion system, which in the case of BCS is a rigorous upper 
bound to the exact many-body value for the given Hamiltonian. Results (b) and (c) 
are also expected to hold for neutral-fermion superfluids (SFs) — such as liquid ^He 
[25, 26], neutron matter and trapped ultra-cold fermion atomic gases [27-38] — where 
the pair- forming two- fermion interaction of course differs from the Cooper/BCS one 
for SCs. 



2. THE COMPLETE BOSON-FERMION MODEL 

The CBFM [12, 13] is described in d dimensions by the Hamiltonian H = Hq-\- 
Hint- The unperturbed Hamiltonian Hq corresponds to a non-Fermi-liquid "normal" 
state, being an ideal (i.e., noninteracting) ternary gas mixture of unpaired fermions 
and both types of CPs namely, 2e and 2h. It is 

ki,si K K 

where as before K = ki + k2 is the CP center-of-mass momentum (CMM) wavevector 
while = fi'^k\/2m are the single-electron, and E±(K) the 2e-/2h-CP phenomeno- 
logical, energies. Here aj^^ («ki,si) are creation (annihilation) operators for fermions 



and similarly (6k) and (ck) for 2e- and 21i-CP bosons, respectively. Two-hole 
CPs are considered distinct and kinematically independent from 2e-CPs. 

The interaction Hamiltonian Hint (simplified by dropping all K ^ terms, as 
is done in BCS theory in the full Hamiltonian but kept in the CBFM in Hq) consists 
of four distinct BF interaction vertices each with two-fermion/one-boson creation 
and/or annihilation operators. The vertices depict how unpaired electrons (subindex 
+) [or holes (subindex — )] combine to form the 2e- (and 2h-) CPs assumed in the 
d-dimensional system of size L, namely 



Hint = L "^/^ ^ f+{k){a+^^.a\^bo + a_kaOk,T^J} 

k 

+L-'^/2 ^/_(A;){a+^a+k,^c+ + a_k,|ak,TCo} (1) 

k 

where k = ^(k-,^ — k2) is the relative wavevector of a CP. The interaction vertex 
form factors f±{k) in (1) are essentially the Fourier transforms of the 2e- and 2h-CP 
intrinsic wavefunctions, respectively, in the relative coordinate of the two fermions. 
In Refs. [12, 13] they are taken as 

[ otherwise. 

One then introduces the quantities Ef and de as new phenomenological dynamical 
energy parameters (in addition to the positive BF vertex coupling parameter /) that 
replace the previous E± (0) parameters, through the definitions 

Ef = ^[E+{0) + E_m and 5e ^ ^[E+{0) - E_m (3) 

where -E'-i-(O) are the (empirically unknown) zero-CMM energies of the 2e- and 2h- 
CPs, respectively. Alternately, one has the two relations 

E±{0) = 2Ef±Se. (4) 

The quantity Ef serves as a convenient energy scale; it is not to be confused with the 
Fermi energy Ep = ^mvp = kBTp where Tp is the Fermi temperature. The Fermi 
energy Ep equals irTi^n/m in 2D and (^^/2m)(37r^n)^/^ in 3D, with n the total 
number-density of charge-carrier electrons, while Ef is the same with n replaced by, 
say, Uf. The quantities Ef and Ep coincide only when perfect 2e/2h-CP symmetry 
holds, i.e. when n = Uf. 

The grand potential for the full H = Hq + Hint is then constructed via 

f](T, L<^, No, Mo) = -keT In [Tre"^^^-'^^^] (5) 

where "Tr" stands for "trace." Following the Bogoliubov recipe [39], one sets 6q , bo 
equal to and Cq , cq equal to ^/Mo in (1), where A^o is the T-dependent number 



of zero-CMM 2e-CPs and Mq the same for 2h-CPs. This allows exact diagonalization, 
through a Bogoliubov transformation, giving [40] 



poo poo 

/ deN{e)[e - /i - E{e)] - 2/csT / deN{e) ln{l + exp[-/?E(e)]} 
Jo Jo 

POO 

+ [E+{0) - 2//]no + ksT / deM{e) ln{l - exp[-P{E+{0) +e- 2//}]} 

Jo+ 

/•OO 

+[2ij- E_{0)]mo + kBT / deM{e)lii{l - exp[-/3{2i^ - E_{0) + e}]} (6) 

Jo+ 

where N{e) and M{e) are respectively the electronic and bosonic density of states, 
E{e) = ^(e-^)2 + A2(e) where A(e) = ^f+{e) + ^/-(e), with no(T) = 
Nq{T)/L'^ and /no(T) = Mq{T)/L'^ being the 2c-CP and 2h-CP number densities, 
respectively, of BE-condensed bosons. Minimizing (6) with respect to A'^o and Mq, 
while simultaneously fixing the total number N of electrons by introducing the elec- 
tron chemical potential namely 

= 0' -^nrr = 0' — = -N (7) 



specifies an equilibrium state of the system with volume and temperature T. Here 
A?" evidently includes both paired and unpaired CP electrons. The diagonalization 
of the CBFM H is exact, unlike with the BCS if, so that the CBFM goes beyond 
mean-field theory. Some algebra then leads [40] to the three coupled integral Eqs. (7)- 
(9) of Ref. [12]. Self-consistent (at worst, numerical) solution of these three coupled 
equations then yields the three thermodynamic variables of the CBFM 

no{T,n,ii), mo{T,n,ii), and l^{T,n). (8) 

Fig.l displays the three BE condensed phases — labeled s+, s— and ss — along with 
the normal phase n, that emerge [13] from the CBFM. 

Vastly more general, the CBFM contains [1] the key equations of all five distinct 
statistical theories as special cases; these range from BCS to BEC theories, which are 
thereby unified by the CBFM. Perfect 2e/2h CP symmetry signifies equal numbers of 
2e- and 2h-CPs, more specifically, nsiT) — niB^T) as well as no{T) — mo{T). With 
(4) this implies that Ef coincides with and the CBFM then reduces to the gap and 
number equations [viz., (11) and (12) below] of the BCS-Bose crossover picture With. 
the Cooper/BCS model interaction — if its parameters V and hujo are identified with 
the BE interaction Hamiltonian Hi^t parameters /^/2fe and fe, respectively. The 
crossover picture for unknowns A(T) and /u(T) is now supplemented by the central 
relation 

A(T) = f^/^) = f^Mif). (9) 

Both A(T) and no(T) and mo{T) are the familiar "half-bell-shaped" order-parameter 
curves. These are zero above a certain critical temperature T^, rising monotonically 
upon cooling (lowering T) to maximum values A(0), ?io(0) and mo(0) at T = 0. 
The energy gap A(T) is the order parameter describing the superconducting (or 




Figure 1. Illustration in the no-mo plane of three CBFfi |oijriiensed phases (the 
pure 2e-CP s+ and 2h-CP s— BE condensate phases and a miyi phase ss) along 
with the normal (ternary BF non-Fermi-liquid) phase n. 

superfluid) condensed state, while no{T) and mo(T) are the BEC order parameters 
depicting the macroscopic occupation that arises below Tc in a BE condensate. This 
A(T) is precisely the BCS energy gap if the boson-fermion coupling / is made to 
correspond to \/2Vhu>D- Note that the BCS and BE TgS are the same. Writing (9) 
for T = 0, and dividing this into (9) gives the much simpler /-independent relation 
involving order parameters normalized in the interaval [0, 1] 

A(T)/A(0) - Vno(T)/no(0) = V^^o(T)/mo(0) ^ 1. (10) 

The first equality, apparently first obtained in Ref. [9] , simply relates the two hereto- 
fore unrelated "half-bell-shaped" order parameters of the BCS and the BEC theories. 
The second equality [12, 13] implies that a BCS condensate is precisely a BE conden- 
sate of equal numbers of 2e- and 2h-CPs. Since (10) is independent of the particular 
two-fermion dynamics of the problem, it can be expected to hold for either SCs and 
SFs. 



3. BCS-BOSE CROSSOVER THEORY 

The crossover theory (defined by two simultaneous equations, the gap and num- 
ber equations) was introduced by many authors beginning in 1967 with Friedel and 



co-authors [41]; for a review see Ref. [42]. The critical temperature is defined by 
A(Tc) = 0, and is to be determined self-consistently with fi{Tc). The two equations 
to be solved, in 2D for the Cooper/BCS model interaction, are [43] 



, f^'^BTc tanhx t; N 1 X , tanha; 

1 = A / ax (it II > nujD)', 1 — X I ax — , (it ji < ncoD) 

Jo ^ -'w^ ^ 

(11) 

* = 1. (12) 



/o exp[e - n{Tc)/2kBTc] + 1 
The last integral can be done analytically, and leaves 

fx{Tc) = A;BTcln(e^^/^^^= - 1). (13) 

The yu(Tc) is then eliminated (numerically) from (11) to give Tc as a function of 
A. Using Tiwd/Ef = 0.05 as a typical value for cuprates, increasing A makes /u(Tc) 
decrease from its weak-coupling (where Tc — > 0) value of Ep down to huiD when 
A ~ 56, an unphysically large value. Fig. 2 displays Tc (in units of Tp) as function 
of A. 



4. GAP EQUATION 

Curiously, the standard procedure in all SC and SF theories of many-fermions is 
to ignore 2h-CPs altogether. Indeed, the BCS gap equation for all T can be derived 
without them. Neglecting in (6) all terms containing mo(T), -E'-(O) and /-(e) leaves 
an 0(T, L*^, |U, A'^o) defining an incomplete BFM. Minimizing it over A'^o (for fixed total 
electron number A^) requires that dQ/dNo = or dQ/dno — 0, which becomes 



^0 



or 



^ ^ 1 + exp{-pE{e)} 

Ef+5e 



2 [E+{0) -2^1] = f jj ' deAr(e)-^ tanh \f3E{^)- (14) 
Using (4) yields precisely the BCS gap equation for all T, namely 

1 = A/ dj ^ — tanh-/jVe^ + A2(r) (15) 

where ^ = e — provided one picks Ef = /x, since A = f'^N{0)/25£ while 6e = Tiwd- 
The companion number equation follows from the last equation of (7) and will thus 
be 

n -n/(T) +2nB(T) (16) 
where nf{T) is the number density of unpaired electrons 

OO 

nf{T) ^ J deN{e)[l - ^ tanh ^/3£;(e)] (17) 
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Figure 2. Critical SC temperatures in units of Tp \or the BCS-Bose crossover 
theory (fuU curve), the BCS value from the exac^mphcit equation (Ref. [21], p. 447) 

1 = A j^^'^ ^'^^^'^'^ dxx~^ in.'oh.x (upper dashed curve) and its weak-couphng solution 
Tc ~ 1.134^a;£) exp(— 1/A) (lower dashed curve). The dot-dashed "appendage" sig- 
nals a breakdown in the BCS/Cooper interaction model when ^(Tg) turns negative, 
as the Fermi surface at jji then washes out. The value of A = 1/2 marked is the 
maximum possible value allowed [45] for this interaction model just short of lattice 
instability. 



while the number density of composite bosons, both with K — {) and with K > 0, is 

oo 

UBiT) = no{T)+nB+{n riB+{T) = j deM{e)^^^^^^^:^—^. (18) 

0+ 

Note that the number equation n = nf{T) + 2nB{T) differs from (12) of the previous 
section which follows from n = Uf (T) only. 

Similarly, ignoring 2e-CPs and keeping only 2h-CPs leads to fl{T, L^, /j,, Mq) 
which to minimize over Mq requires that dfl/dMQ = dQ/dmo = and, since E{^) = 
E{—^). This again leads to (15) but now with the companion number equation 

m^nf{T) -2mB{T) (19) 

with the same previous nf(T) and where 

oo 

TUBiT) = mo{T)+mB+{n ruB+iT) = J deM{e) _ ^ . (20) 

0+ 



However, ignoring either 2e- or 2h-CPs does not give the entire BCS ground-state 
energy, as we now show. 



5. CONDENSATION ENERGY 



The T = condensation energy per unit volume according to the CBFM (i.e., 
with both 2e- and 2h-CPs) is 

Es-Er, _ ns{T = o)-nn{T = o) 

Ld - Ld ^^^^ 

since for any T the Helmholtz free energy F = Q + = E — TS, with 5" the 
entropy, and /x is the same for either superconducting s or normal n phases with 
internal energies Eg and E^, respectively. In the normal phase tt-q = 0, mo = so 
that A(T) = for all T > 0, and (6) reduces to 



»n(r = o) 



poo pfi pO 

/ deN{e)[e-f,-\e-f,\] = 2 / deN{e)[e- f,] = 2 / d^N{0^. (22) 

Jo Jo J-fj, 

For the superconducting phase at T = 0, and when no{T) — mo{T) and nsiT) = 
mB{T) hold, one deduces from (4) and (6) that /i = Ef. Letting A(T = 0) = A in 
(6) and putting Se = tkoo while using (4) gives for the superconducting phase 



= 2hwDno{Q) + 2 f d^NiO^ - 2 / "° d^N{0./eT^. (23) 
J-n Jo 

Subtracting (22) from (23) and putting A^(^) = -^(0), the density of electronic 
states at the Fermi surface, one is left with 



T71 T71 rfi'jJ 



d^ - v^eM^j (CBFM). (24) 



Ld 

Employing Eq. (2), p. 158 of Ref . [44] the integral becomes 



2 2 ^ ^ 2 nu;D + y(^i:^^FTA2 



W^n(^)-W-1^^0l^]. (25) 



A^o 2 \2nujD J 4 16(na;o)2 \[huD 



Using (9) for T = and weak coupling / — > implies that A = />/no(0) = 
/\/too(0) ^ so that (24) yields the expansion 



Eg — En 



Ld— ^ 2hujDno{0) 



+ 2Ar(0) 



- A^ In — A^ — + O 



2 "^\2hujDj 4 16{nuD)^ [[huot 



(CBFM). (26) 



Given that for small A 



sinli(l/A) A^o 



2TujJd exp(— 1/A) 



the log term in (26) is just 



/ A \ ^ 2?iWD 



since 



A = VN{0) 



so that (26) finally simplifies to 



Es — En 



-iV(0)A2 



1 / A 

1+ T 



O 



D 



(27) 



(28) 



(29) 



(CBFM). (30) 



By contrast, the original BCS expression from Eq. (2.42) of Ref. [6] is 

Eg — En 



N{0){nu}Dy 



l-Jl + iA/houDY 



(BCS) 



(31) 



which on expansion leaves 



Eg — En. 



iiV(0)A^ 



^ 1 

4 XhUD 



+ 



huj 



D 



(BCS). (32) 



Thus, the CBFM condensation energy (30), and consequently its ground-state en- 
ergy, is lower (or larger in magnitude) than the BCS result (32). Therefore, the 
CBFM satisfies a prime expectation of any theory that improves upon BCS, which 
being based on a trial wave function gives a ground-state energy which is a rigorous 
upper bound to the exact energy associated with the BCS Hamiltonian ground state. 
Consequently, there is no a priori reason why the CBFM is limited to weak coupling, 
at least for all A < 1/2 [45]. 

What happens on ignoring either 2e- or 2h-CPs, as seems to be common practice 
in theories of SCs and SFs? Starting from (6) for T = 0, and following a similar 
procedure to arrive at (23) but without 2h-CPs such that /_ = 0, mo(0) = and 
no(0) = AV/2, one gets 



/O rhuJn , . 



Subtracting (22) from (33) gives 



Eg — En 



Ld 



/•huJn 

= hojonoio) + iv(o) J (e - Ve' + a^ j 



(33) 



(34) 



which is just one half the Ml CBFM result (24). Furthermore, if [{Eg — En)/L'^]- 
is the contribution from 2h-CPs alone, assuming now that /_|_ = and no(0) = we 
eventually arrive at precisely rhs of (34) but with mo(0) = A^//^ in place of no(0). 
Hence 



Eg — Er. 



Eg — Er. 



1 



1 + 



1 

4 \huj 



D 



+ 



D 



(35) 



which again is just one half the full CBFM condensation energy (30) that in leading 
order in A was found to be the full BCS condensation energy. 

Including both 2e- and 2h-CPs gave similarly striking conclusions on generalizing 
[22, 23] the ordinary [18] CP problem from unrealistic infinite-lifetime pairs to the 
physically expected finite-lifetime ones. 



6. CONCLUSIONS 

The recent "complete boson-fermion model" (CBFM) contains as a special case 
the BCS-Bose crossover theory which, at least for the Cooper/BCS model interaction, 
predicts virtually the same TcS to well beyond physically unreasonable values of 
coupling than the allegedly less general BCS theory where the number equation is 
replaced by the assumption that fx = Ep ■ 

The CBFM reveals that, while the BCS gap equation for all temperatures follows 
rigorously without either electron or hole pairs, the resulting T = condensation 
energy is only one half the entire BCS value. In view of this, if BEC is at all 
relevant in SCs and SFs taken as many-fermion systems where pairing into bosons 
undoubtedly occurs, two-hole CPs cannot and must not be ignored. 
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